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ABSTRACT 

We present the linear inverse redshift space operator which maps the galaxy density 
field derived from redshift surveys from redshift space to real space. Expressions are 
presented for observers in both the CMBR and Local Group rest frames. We show 
how these results can be generalised to flux-limited galaxy redshift surveys. These re- 
sults allow the straightforward reconstruction of real space density and velocity fields 
without resort to iterative or numerically intensive inverse methods. As a corollary to 
the inversion of the density in the Local Group rest frame we present an expression 
for estimating the real space velocity dipole from redshift space, allowing one to es- 
timate the Local Group dipole without full reconstruction of the redshift survey. We 
test these results on some simple models and find the reconstruction is very accurate. 
A new spherical harmonic representation of the redshift distortion and its inverse is 
developed, which simplifies the reconstruction and allows analytic calculation of the 
properties of the reconstructed redshift survey. We use this representation to analyse 
the uncertainties in the reconstruction of the density and velocity fields from red- 
shift space, due to only a finite volume being available. Both sampling and shot-noise 
variance terms are derived and we discuss the limits of reconstruction analysis. We 
compare the reconstructed velocity field with the true velocity field and show that 
reconstruction in the Local Group rest frame is preferable, since this eliminates the 
major source of uncertainty from the dipole mode. These results can be used to trans- 
form redshift surveys to real space and may be used as part of a full likelihood analysis 
to extract cosmological parameters. 
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1 INTRODUCTION 

Galaxy redshift surveys are distorted by the intrinsic mo- 
tions of galaxies along the line of sight. This distortion leads 
to a number of effects. On large, linear scales the radial dis- 
placement of galaxy positions is enhanced, leading to com- 
pressed structures along the line of sight. On smaller scales, 
near turnaround, the structure will appear completely col- 
lapsed along the line of sight, and thereafter will begin to 
invert itself. Finally on very small scales, after virialisation, 
the random motions of galaxies in clusters will stretch out 
the structure in redshift space. 

These distortions can complicate measurements of the 
statistical properties of the large scale galaxy distribution, 
and can lead of confusion in the interpretation of redshift 
maps. However, the distortions also contain valuable infor- 
mation on the structure of cosmological velocity fields. Fur- 
thermore, in gravitational instability theory the parameter 
controlling the degree of distortion is a function of the mean 



density of the Universe (see the excellent review by Hamil- 
ton, f997, on linear redshift distortions). 

One of main problems of cosmography is how to map 
the large scale structure seen in galaxy catalogues from red- 
shift space to real space. In recent years there has been much 
interest in reconstructing redshift surveys, with the advent 
of the IRAS all-sky surveys, QDOT and 1.2Jy, and in partic- 
ular the new Point Source Catalogue redshift survey (here- 
after PSCz; see Saunders et al 1996). While it is straightfor- 
ward to write down the transformation from real space to 
redshift space galaxy distributions, the inverse relation has 
proven more elusive. 

The standard approach to reconstructing real space and 
the velocity fields has been to apply some iterative method, 
to solve for the density field (Kaiser et al 1991, Taylor & 
Rowan-Robinson 1993, Branchini et al 1998) or reconstruct 
galaxy positions, again iteratively (eg, Yahil et al 1991). Pee- 
bles (1989) showed that the trajectories of galaxies could be 
solved numerically by a Least Action principle and Shaya 
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et al (1993) applied this to large scale structure. Recently 
Croft & Gaztanaga (1997) have shown that applying the 
Zel'dovich approximation to a Least Action principle leads 
to a combinatorics problem of finding the least path between 
initial and final conditions. Valentine et al (1999) have gener- 
alised this to selection function weighted galaxy surveys and 
applied the method to the Point Source Redshift Catalogue. 
Again, the process is iterative and may be time consuming. 
Going to a higher order reconstruction appears to be a slow 
and painful process. 

An alternative method has been to expand the density 
field in spherical harmonics and solve for each multipole, ei- 
ther from a numerically intensive inversion method (Fisher 
et al 1995, Zaroubi et al 1995, Webster et al 1997), or by 
numerically solving the redshift space Poisson equation to 
find the velocity potential field (Nusser & Davis 1994). The 
former method has the advantage of calculating the uncer- 
tainties in the reconstruction, but assumes a Gaussian dis- 
tributed density field and requires regularization in the form 
of a Wiener filter. The latter has the advantage of being 
extendable into the Zel'dovich regime, but again must be 
solved numerically if flux-limited. 

Finally, Tegmark & Bromley (1995) proposed a Green's 
function solution to the linear inversion problem, allowing 
one to regain the real space density and velocity fields from 
a straightforward convolution. However their method was 
only valid for volume limited galaxy distributions. 

With the advent of the PSCz and forthcoming comple- 
tion of the AAT 2dF (Maddox et al 1998) and Sloan Digital 
Sky Survey (Szalay et al 1998) an accurate, but straight- 
forward method for removing the distortion in flux-limited 
redshift surveys is required. Reconstruction can be seen as 
required for comparison of velocity field data with the in- 
ferred flow in redshift surveys, or a general method for cor- 
recting the distortion before subsequent analysis. 

Since redshifts are taken in the rest frame of the ob- 
server, the inversion method should deal with the local 
dipole motion of the observer. The dipole itself is of in- 
terest, since a reconstruction of the real space dipole can 
be compared with the observed CMB dipole and again put 
constraints on the mean density of the Universe (Rowan- 
Robinson et al 1990, Kaiser & Lahav 1988, Rowan- Robinson 
et al 1999, Schmoldt et al 1999). Since the dipole is nonlocal, 
the whole volume of the redshift survey usually has to be 
reconstructed in order to estimate just the dipole. 

In this paper we address the problem of linear recon- 
struction of galaxy redshift surveys, and put it on a firm 
theoretical foundation. In Section U we derive the inverse 
spherical redshift distortion operator by considering the dis- 
placement fields of galaxies in real and redshift space. This 
operator allows one to transform directly between the lin- 
ear redshift density field and the linear real space density 
field, without iteration. This allows a quick and simple, yet 
accurate, method for removing the linear distortion seen in 
redshift maps. We derive versions of this inverse operator for 
volume and flux-limited redshift surveys and for observers in 
the Cosmic Microwave Background (CMB) and Local Group 
rest frames. As a corollary of the Local Group inverse red- 
shift operator we derive an expression for estimating the real 
space observer's dipole motion directly from a redshift sur- 
vey, without a full inversion. This removes the effects of the 



redshift distortion and corrects for the well-known "Rocket 
effect" (Kaiser & Lahav 1988). 

We develop a spherical harmonic representation of the 
distortion operator and its inverse in Section [| This allows 
both practical and convenient application of the operators to 
redshift surveys, and simplifies the calculation of the prop- 
erties of the reconstructed density and velocity fields. We 
demonstrate the method by reconstructing both linear den- 
sity fields and the real space dipole from a simple model of 
a redshift survey. 

In Section |j| we use our spherical harmonic formalism to 
estimate the uncertainties in the reconstruction of the den- 
sity and velocity fields with only a finite survey volume. We 
derive the statistical properties of the reconstructed veloc- 
ity field, taking into account cosmic variance and shot noise. 
This allows us to place limits on the accuracy to which one 
can do such reconstructions. 

We summarise our results and present our conclusions 
in Section |H| 



2 MAPPING FROM REDSHIFT— SPACE TO 
REAL-SPACE 

The problem of finding an inverse redshift-space operator 
depends on the specifics of the redshift survey to which it is 
to be applied. The redshift survey may be volume- or flux- 
limited. In the case of the latter, the mean observed density 
of galaxies is a function of distance from the observer and 
requires an extra correction. The redshifts of the galaxies 
in a survey are taken in the observer's rest frame, and so 
include the motion of the Earth with respect to the CMB. 
Traditionally this has been transformed to either the Local 
Group rest frame - where the Local Group represents the 
largest gravitationally bound structure to which we belong, 
and whose size is close to the linear regime - or the CMB 
rest frame. In the rest of this paper we shall refer to the 
"Local Group" rest frame in rather a loose way, referring to 
a hypothetical observer moving in the local linear flow field 
with respect to the CMB rest frame. 

In the next few sections we find solutions to the inverse 
dist ortion in the CMB frame for volume-limited (Section 
2.1) and flux-limited (Se ctio n 2.3) surveys, and in the Lo- 
cal Group frame (Section 2.4). In addition the effects of the 
"Local Group" motion on redshift space and its inverse op- 
erator lead us to an expression for the calc ulati ng the real 
space dipole from a redshift survey (Section 2.4). 

We begin by considering the linear displacement of 
galaxies and derive the inverse redshift operator for a volume 
limited redshift survey. 



2.1 Volume limited redshift inversion 

The forward problem, of mapping a density field from real 
space to redshift space, can be solved using the transforma- 
tion^] 

s = r + (r.v)r (1) 



* Throughout we shall use velocity units for distances. This 
is equivalent to setting the Hubble parameter Ho = 1 in the 
formulae. 
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where s and r are the comoving redshift and real-space co- 
ordinates, and v is the peculiar velocity field. In linear the- 
ory the velocity field is calculated from the density field by 
(Peebles 1980) 



v(r) 



/(n„ 



4tt 



d¥ <5(r)- 



(2) 



where /(fl m ) — f^m 6 is the growth rate of the density field 
and Q m is the mass density parameter of the Universe. 

It is useful to introduce a displacement vector field, £, 
defined by 

dt V 

which denotes the displacement in real space from some ini- 
tial, or Lagrangian, position to the present day. In linear 
theory the velocity is proportional to the displacement, 



(3) 



V = /(Om)£. 



(4) 



To allow for the possibility of linear biasing of the galaxy 
density field we introduce the linear bias factor, b, defined 
by S g = b5 m . If velocities are estimated from the galaxy 
distribution, rather than the true density, the constant of 
proportionality between velocities and displacements is j3 
where /3 = 

The real space and redshift space coordinates, r and s 
respectively, can be related to a set of initial or Lagrangian 
coordinates, q, by 



r = q + £, 
and 

a = q + f. 



(5) 



(6) 



The redshift space displacement field, £ s , is related to the 
real space displacement, £, by (Taylor & Hamilton 1996) 



& = Piiti 
where 



(8) 



is the redshift space projection tensor, and 6^ is the Kro- 
necker tensor. From continuity we can find the linear redshift 
density by taking the divergence of the redshift displacement 
field: 

5 s = -V.| s . (9) 
Substituting equation (Q) into equation (^) yields 

5 s = -(V.£ + /?(f i f j V i & + 2r- 1 f i &)) 

= [1 + /3(d 2 +2r- 1 d r )V 2 ]S (10) 

where we have used the real space continuity equation 

S = -V.£, (11) 

and V -2 is the inverse Laplacian. To linear order spatial 
derivatives in redshift space are equal to derivatives in real 
space, and the coordinates of fields are the same in both real 
and redshift space, i.e. 8(s) — S(r). 

Equation ( |lo| ) can be written much more succinctly if 
we introduce the linear, spherical redshift space operator 
(Hamilton 1997) 

S = l + /3(a 2 + 2r~ 1 9 r )V" 2 , (12) 



and write the transformation from real to redshift density 
fields as 



5 s = SS. 



(13) 



Having formulated the redshift distortion in the lan- 
guage of displacements the inverse operator, is now 
straightforward to find. We begin by inverting equation (pj), 
relating the real and redshift displacement fields, 



where 



1 + /3 



(14) 



(15) 



is the inverse redshift projection tensor. Here we have used 
the identity s = r. Taking the divergence of both sides of 
equation ( |l4| ) will give us the inverse relationship between 
redshift and real space density perturbations, 



: (d 2 s +2s- 1 ds)\7- 2 ]S s 



(16) 



(17) 



1+/T 

The inverse spherical redshift operator is given by 
S- 1 = 1 - Y^p( d ° + 2s -1 &)V- 2 

To invert the redshift distortion we have used only the con- 
tinuity equation in real and redshift space, and assumed po- 
tential flow in bo th t he real and the redshift displacement 
fields. In Section 



2.2 



we shall show that if the linear real 
space displacement field is curl-free then so is the linear 
redshift displacement field. 

Having found the inverse operator we now verify that 
it satisfies the relation 



SS' 1 = 5 _1 S' = 1. 



(18) 



In order to show that equation ( |17| ) does satisfy this it is 
useful to first find the relation between S and Vij, and S~ 
and 'P~ j 1 . From the redshift continuity equation we find 



-V.| s 



(19) 



Expanding 6 in the same way we find that S and V and S 1 
and V~ x are related by 

S = (V^V^V" 2 , S- 1 = (V^V^V- 2 , (20) 

Multiplying these operators together yields 

SS' 1 = (V j PyV i )V- 2 (V fc 7' fe l 1 Vi)V- 2 , 

= V l P lJ (V J V- 2 V fc )m 1 V ; V- 2 (21) 

where we have highlighted the operator VjV _2 Vfc in the 
second line which is applied to the real space displacement 
field, £. If we apply this operator to a general potential field 
Ai = V%4> we find 



V,:V~ 



ViV" 
Ai. 



(22) 



Hence for a potential vector field this operator is the identity 
matrix 



(23) 



Substituting this into equation ( pl[ ) completes the proof 
SS' 1 = ViVijiVjV^V^V^ViV- 2 
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VAfv ; v~ 2 
1. 



(24) 



2.2 Is the linear redshift displacement field 
irrotational ? 

To establish the same proof for S~ 1 S we need to show that 
the operator ViV~ 2 Vj, which this time is applied to £ s , is 
acting on a potential field. Introducing the redshift vorticity 
vector, u> a defined as 

u s = /3V x (25) 
and given that the real space displacement field is curl-free, 
w ee V x £ = 0, (26) 
we find 

= /3V x f(f .£) 

= /?(f.£)(V x f) + /?(f x V)(f.£) 

= 0. (27) 

Hence the linear redshift displacement field is curl-free and 
our proof holds. 

2.3 Flux limited redshift inversion 

Since most redshift surveys are flux limited rather than vol- 
ume limited it is useful to find the flux limited linear redshift 
operators. The density of galaxies seen in a redshift survey 
can be represented in real space by 

p = 4>(l + S), (28) 

while in redshift space the galaxy density can be written as 

p s =4> S (1 + S S ). (29) 

The mean observed density of galaxies in the survey is given 
by the selection function, <j), in real space and s = <fi(s) in 
redshift space. Conservation of galaxy numbers imply that 

p(r) d\ = p s (s) d A s (30) 

and so the real and redshift density fields are related by 

p = p s det£V, (31) 

where 

dsi 



drj 



(32) 



is the distortion tensor mapping real to redshift space. Ex- 
panding the determinant of the distortion tensor to first or- 
der we find 



dct Vi 



1 + 



1+, 



where we have used the relationship 



: n + 



p 



(33) 



(34) 



to express the expansion in terms of the redshift displace- 
ment vector. We can also expand the real space selection 
function in terms of redshift space quantities: 



■(s4 3 )ds ln< 



(35) 



1+13' 

Combining these expressions we find to first order that the 
spherical redshift operator for a flux limited galaxy survey 
is 



5=[l-JL-(d 2 s +a'(s)s- 1 d s )V- 2 ]5 s 



where 



q'(s) 



dlns 2 (f>- 1 (s) 



(36) 



(37) 



dins 

is related to the local slope of the selection function. Hence 
the inverse redshift operator for flux limited redshift surveys 
is 



s-' = i 



P 



1 + , 



; (d 2 s +a'(s)s- 1 d s )V- 



Similarly the forward redshift operator is 
S = l+P{d 2 + a{r)r- 1 d r )V- 2 
where 

. , d\nr 2 <h(r) 
air) = — — — 
dlnr 

a and a' are related by 
a(r) + a'(r) = 4. 



(38) 



(39) 



(40) 



(41) 



It is clear from these relations that the effect of changing 
from a volume limited redshift survey to a flux-limited red- 
shift survey can be easily incorporated by the substitution 
2 — > a or 2 — > a' in the forward and inverse operators respec- 
tively. This simple transformation will allow us to generalise 
results from volume- to flux-limited. 



2.4 Inverse Redshift Operator in the Local Group 
Frame 

So far we have only considered the distortion operator and 
its inverse from the point of view of an observer at rest 
with respect to the CMB. However the true rest frame of an 
observer is influenced by local, nonlinear gravitational in- 
teractions. Since nonlinear motions induced by local galax- 
ies cannot be treated in the linear framework it is usual to 
transform to the linear part of our motion, that of the Lo- 
cal Group of galaxies (whether our Local Group also has a 
substantial nonlinear component to its motion is an issue we 
shall side-step here). One can then either work directly in 
the Local Group frame, or use the known dipole, measured 
by the CMB dipole anisotropy and transform to the CMB 
frame. However, there are a number of reasons for work- 
ing in the Local Group rest frame. Firstly the linear theory 
formalism is only valid so long as the perturbative variable 
r.(v(r) - v(0))/r when r — » 0. In the CMB frame the 
velocity field does not necessarily vanish at the origin, and 
an extra, un-necessary assumption is needed to satisfy this. 
In addition, if our relative motion does not vanish at the 
origin, divergences appear in the analysis of our dipole, cal- 
culated from redshift surveys. We discuss this effect in detail 
below. Finally, as we shall show in Section 3 reconstruction 
of the velocity field in the Local Group frame reduces the 
uncertainties due to finite survey volumes and shot-noise. 
The transform from real space to redshift space in the 
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Local Group rest frame can be found by subtracting the 
Local Group displacement from the redshift displacement 
vector: 



s,LG 



Viiti - m&M = - Va ~ *«)&(°)- (42) 



Super- or sub-script LG denotes a variable measured in the 
Local Group frame, and £(0) is the Local Group displace- 
ment, or dipole. This is equivalent to using r.(v — v(0)) in- 
stead of r.v as the redshift coordinate. From equation ( |42| ) 
the transformation to redshift space is straightforward. Tak- 
ing the divergence of the redshift displacement field in the 
Local Group frame we find 



[1 + P{d 2 r + 2r-\d r - oWo))V- 2 ]<5 



(43) 



where 9 r= o is the radial differential operator at the origin 
and <9 a =oV -2 = — J d 3 r {s .r) / Anr 2 is the radial velocity op- 
erator at the origin. It is useful to note that ViSiSj£j (0) = 
£j(0)(ViSjS;) since £(0) is a fixed vector at the origin. 

In the Local Group frame the redshift operator is then 



i + /3(a 2 + 2r- 1 (a r -a r=0 ))v- 



(44) 



Finding the inverse is slightly trickier. A straightforward in- 
version of equation (W2) yields 



6 = V^' iG + (*«-^«%(0) 



(45) 



Therefore to find the real space displacement field we need 
to know what the Local Group displacement is. One could 
take the known dipole from the CMB here, but to calculate 
the effect self-consistently the Local Group motion should 
be estimated from the redshift catalogue itself. 

One might suppose that equation (0) would help us 
here, since it transforms real displacements to redshift dis- 
placements. But if we let r -> we find £ S ' LG (0) -> £(0). 
This is just a statement that the only displacement not af- 
fected by the redshift transformation is the observer's, since 
this is by definition the origin of coordinates in both systems. 
But while we have demonstrated that the redshift and real 
space Local Group displacements are equal, this does not 
imply that the Local Group motion estimated from redshift 
catalogs are the same in real and redshift space, and it is 
the latter quantity we need. To illustrate this, consider the 
dipole generated at the origin. This can be calculated from 
the density field by 



(46) 



Similarly, one can calculate the redshift dipole (Kaiser & 
Lahav 1988) 



rco) 



d s r r 
An r 2 



5» 



d r r 

An r 2 



+p I ^-^{dl + 2r- 1 d r )v- 2 s, 



(47) 



which is clearly different from the true dipole. If we only 
consider the first term, which arises from the intrinsic dipole 
motion of the observer, 



€*(o)=#(o)| If, 



(48) 



we find a logarithmic contribution to the dipole from the 
survey geometry. This is the well known "Rocket Effect" 
(Kaiser & Lahav 1988). This contribution is generated by 
the reflex action of the distorted survey on the measured 
dipole. The divergence of the estimated dipole arises when 
points in real space are mapped to the origin in redshift 
space. In practice we would not expect this divergence to 
appear if the velocity field is coherent over some scale. In 
this case the local matter field will be moving at the same 
velocity as the observer and will not be mapped to the ori- 
gin. Hence the way to avoid such divergences is to smooth 
the velocity field on some scale. In practice we will have to 
smooth on large scales to allow the use of linear theory (see 
Section 3.6 below). 

A second divergence appears at large radii. As r — » 
oo the Rocket term diverges logarithmically, so the redshift 
dipole in the Local Group frame does not converge (Kaiser 
& Lahav 1988). This consideration shows us that the true 
dipole can only be estimated over a finite range of radii. In 
practice a reconstruction should be done for different survey 
depths to test if the true dipole has converged. 

The true dip ole can be estimated from redshift space 
via equation ( |46p and the divergence of equation ( fiB] ) 



m = 



-j- 4 V -£M 

4-7T r z 



da s 



|3. s 
47T S 2 

2/? 
3(1+ /J) 



/? 



VP 
dr 
r 



{d 2 + 2s~ 1 d s )V- 



(49) 



Solving this we arrive at the dipole solution 







1+0 



where 



A- 



23 



3(1 + /?) 



dr 
r 



(d 2 + 2s- 1 d s )\7~ 2 ]Sl G , (50) 



(51) 



The presence of a logarithmic term, J dr/r, may seem wor- 
rying but it is there to cancel the divergent term that arises 
from estimating the dipole in redshift space from the Lo- 
cal Group frame. Again the divergence at small scale can 
be removed by smoothing the velocity field, or estimating 
the dipole for a ball of matter at the origin. The redshift 
density field in equation ( |5o| ) is acted on by the inverse op- 
erator, but no correction is made for the dipole, so this field 
is not a mapping to the real space density. We present the 



correct expression for this below. In addition, in Section 3.4 
we show that the dipole is coupled only to the redshift space 
mass dipole. Hence we can estimate the Local Group dipole 
directly from redshift surveys, without reconstructing the 
whole survey. 

The inverse redshift distortion operator in the local 
group frame is somewhat more complicated than in the 
CMB frame, and can be written as 



f=[l- 



1 + /3 



(d 2 + 2s- 1 d a )v- 2 }s s LG 



3 j-s.m, (52) 



1+/3 s 
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where £(0) is estimated from equation (|5Cj), or taken from 
the CMB dipole. 

Writing this in the form of an operator we find that the 
inverse redshift operator in the Local Group frame is 



a-l 



1 + — ^-T"»=0V 



1 + /3 As 



x 1 







1+0 



(a s 2 + 2s" 1 9 s )V" 



(53) 



These operators are valid for volume-limited redshift 
surveys. For flux limited surveys the effects of the selection 
function are easily incorporated by making the substitution 
2 — > a(r) in the operator S and substituting 2 — > a'(s) in 
the operator S -1 , and noting that these are now functions 
of radius and so should appear to the right of the integral. 
This completes our task of finding the linear inverse redshift 
distortion operators. 



3 SPHERICAL HARMONIC 

REPRESENTATION OF THE REDSHIFT 
OPERATORS 

3.1 The redshift space operators 

In the case of spherical survey geometry and radial distor- 
tions the most convenient basis is that of spherical harmon- 
ics (Heavens & Taylor 1995, Ballinger at al 1995, Fisher at al 
1995, Tadros et al 1999). In this basis we can use the prop- 
erties of the eigenvalues of the Laplacian to transform from 
an integro-differential operator to an algebraic operator. We 
shall assume for convenience that the survey is all sky. For 
a spherical harmonic function, tpe m (k), defined by 



d 3 rip(r)ji(kr)Yl m (r) 



the Laplacian becomes 



V 2 V> = ( dr + ~d r - ^ ) i> = -k 2 il> 



(54) 



(55) 



where L 2 = £(£ + 1) is the square of the angular momentum 
operator. In this basis we can rewrite the distortion operator 
as 



s=i+p[ 1 



L 

k 2 r 2 



whose inverse is 



s -x 



1+0 



(56) 



(57) 



These operators contain both wavemodes and radial coordi- 
nates and should be used in expressions such as 

S ( r ) = \/f / dkk 2 J2 5 lm(k)S- l (k,r)ji(kr)Y e * m (r), (58) 

transforming real space harmonic modes into the redshift 
space density field. The physical interpretation of these op- 
erators is straightforward. For S an isotropic term propor- 
tional to is added to the true field and the angular part 
proportional to L 2 is subtracted, leaving only a radial distor- 
tion. The inverse operator, S' 1 , replaces the angular term, 



producing an isotropic distortion, and then reduces the field 
by a factor 1/(1 + p). 

To including the effects of a selection function we add 
an extra term 



1+0 1 



k 2 r 



L \ a- 2 2 

+ OrV , 



(59) 



where the integro-differential term can be replaced using the 
recurrence relation for spherical Bessel functions: 



-<9rV ji{kr)Yt m (f) 



{ij^jkr) -{£ + !) j l+1 {kr)) 
(2£ + l)kr 



Y tm {r) 



causing some mixing between ^-modes. Moving to the Local 
Group frame requires another extra term 



l+0\ 1 



L 

k 2 r 2 



+ 







'Sirkr 



&nSmO I d 3 r' S D (r') (61) 



where the last term is the dipole contribution in spherical 
harmonics (see section 4.5), and the integral operator acts 
to shift positions to the origin. 



3.2 The transformation properties of harmonic 
modes 

Another representation is the transformation of the spheri- 
cal harmonic modes of the redshift survey. In this case the 
modes transform according to the relation 

S? m (k) = dk'S tm {k')Kt{k',k)(62) 

where we have used the integral relation (Watson 1966) 



drjt(kr)jt(k'r) = ^ - ^ K t (k', k) 
and the radial function Ki(k! , k) is defined as 



(63) 



z > z 



z > z. 



(64) 



The inverse transformation, from <5| m (fe) to 8t m (k), is done 
in the same way, with the substitutions (1 + 0) — > (1 + /3) _1 
in the first term in equation (p2[), and /3 — > —(3/(1 + (3) in 
the second. 

Similarly the redshift transformation in the Local 
Group frame is 

s tt G i k ) = (1 + 0)hm{k) - ^jj- f dk' 8 lm (k')Kt(k', k) 



+lp6& I dk' S em (k')(k'/k 2 ) 



(65) 



where the radial transform of the dipole term is completed 
by use of the integral (Watson 1966) 



nfrfcffaQ- v^ r[(l + 2)/2] 
rdr 3l {kr)- fc2 r[(£ + 1)/2] 



(66) 



In Figure |l| we show set of Gaussian density fields, 
5t,m{r) in real and redshift space, using the transformation 
equation (|62[). The harmonic modes of a constrained Gaus- 
sian field can be generated by the relation 
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Figure 1. £ = 1 to 6 and m = density wave in real space (black 
line), redshift space (light-grey line) and recovered density field 
(dark— grey, dotted line) in the CMB frame with {3 = 1. 



5im(k) = ^P{k)e w ^ k) 



(67) 



where P(k) is the linear power spectrum derived by Pea- 
cock & Dodds (1994) and 8i m (k) is randomly distributed 
for each mode. The waves have I = 1 to 6 and m = 0. 
Since the angular modes are independent and the distortion 
is m-independent this choice of azimuthal mode does not 
affect this demonstration. The real space density wave was 
recovered using equation ( p8| ) . It is clear in this test of lin- 
ear reconstruction that modes can be recovered with great 
accuracy. 



3.3 The redshift power spectrum 

Having calculated the harmonic modes of the density field in 
redshift space, we can calculate their correlations. This has 
been calculated before by Heavens & Taylor (1995) for a dis- 
crete set of spherical harmonics and by Zaroubi & Hoffmann 
(1996) for a set of Fourier harmonics. However the methods 
presented here provide a simpler expression. 

Taking the two-point expectation value of the harmonic 



modes we find 

<«m(*)#m(*')> = {l+P) 2 P(.k)k- 2 5 D {k-k') 



(2£ + l) 2 J 



dkiP(ki)k 1 Ke(ki, k)Ke{ki, k') 



-0(1 + l3) £ -^Y(P(k)k- 2 K t (k', k) + P{k')k'-\ e (k, k% 

(68) 

and zero for different £'s or ms. As well as redshift space dis- 
tortions creating correlations between different modes, the 
second term shows that the shape of the redshift power spec- 
trum depends on the overall shape of the real power spec- 
trum through a convolution. 



3.4 The cosmological dipole 

In Section 2.4 we derived an expression for the real space Lo- 
cal Group dipole. Again this takes a simple form in spherical 
harmonics. Expanding 5 s in harmonics, and using the rela- 
tion 



-Y u (r) 



(69) 



and the orthogonality of spherical harmonics, we can reduce 
equation (pO) to 



dkk 2 drSj^ifar^tiityjiikr) (70) 



6(0) 



The radial integral can be performed to give 
1 



6(0) = / dkkSUk) 

V6ttM(1 + 13) J 

x\j (t) + 2/3/3(j (t)+h(t)/t- 



«(*))]& (71) 



where R and ro are the upper and lower radial limits of the 
survey and Ci(z) = — J z dtcos(t)/t is the cosine integral. 

Alternatively one can express the dipole in terms of the 
redshift space density field without reference to the har- 
monic expansion. Integrating equation (^5|) and using the 
definition of the density dipole we find the real space dipole 
reduces to 

j., s 1 fd 3 s «/ n- ( 1 0/1 1 2 

^ } = mxW) J 3 ( 3 } 3 b + 3 b + & — s 

Equation ([n]) or (^2|) can also be used to estimate the real 
space dipole contribution in shells. In the limit of no dis- 
tortions (f3 — > 0) equation (j?^) reduces to the real space 
dipole, given by equation (^), and equation ( |7l| ) reduces to 
the harm onic equation for the dipole (equation (109), Sec- 
tion 4.5). 



(72) 



Hence we see that the real space velocity dipole can be 
reconstructed from the redshift space density field, without 
a full reconstruction. Added complications arise if there is 
incomplete sky coverage, since modes will be mixed, and this 
mixing must be included in a dipole estimate. 

Finally, we can add on the effects of a selection func- 
tion. Making the substitutions 2 — *■ a in equation ( |5o| ) and 
working through, wc find the selection function correction 
term to equation ([72]) is 



it (o) = 



d s S s (s)s 



^2 9 " ln< 



(/). (73) 
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Figure 2. True CMB redshift and recovered dipole from a ran- 
dom Gaussian realisation of a density field, calculated in the CMB 
rest frame. Only the amplitude of the dipole is plotted. The light 
grey line is the redshift space dipole. The darker lines are the 
true and reconstructed dipoles. The fields have been Gaussian 
smoothed on a scale of 5/i -1 Mpc. 



where the factor of 2 in A is also transformed to a' . 



3.5 Simulated reconstruction of the dipole motion 

Figure |2| shows the amplitude of the dipole reconstructed 
from a simple random Gaussian model of the density field. 
The true dipole was calculated from equation (fn]) with (3 = 
(black line) . Transforming the spherical harmonic modes to 
redshift space in the CMB frame was then done via equation 
(^), with /3 = 1, and the dipole re-calculated in redshift 
space (grey line). Finally the true dipole (dark grey line) 
was recovered via equation (fnj). It is clear from figure g 
that the true dipole can also be recovered to good accuracy 
from equation (fri]). 



clusters can also lead to divergences in the calculated veloc- 
ity field when the distance between galaxies becomes small. 
Hence smoothing would seem to be a good thing. 

However smoothing in redshift space is not the same as 
smoothing in real space, since the coordinate system changes 
between the two. An isotropic smoothing kernel in one space 
will not be isotropic in the other. Anisotropic smoothing 
kernels, such as in adaptive, or Lagrangian smoothing which 
uses the local moment of interia as a kernel, will preserve the 
shape of structure and can be mapped back to real space 
smoothed with the correct kernel. 

The spherical harmonic decomposition uses sharp fc- 
and (£, m)-space filtering. This has the advantage that both 
redshift and real space will have the same smoothing, al- 
though since the harmonics are related by a convolution over 
fc-space we need to sample higher fc-modes in the inverted 
space. This will be limited by nonlinearities in the higher fc 
range. The sharpness of the fc-space filter can be removed 
by introducing extra filters. 



4 PROPERTIES OF RECONSTRUCTED 
REDSHIFT SURVEYS 

Having found the inverse redshift space operator, we can 
now transform redshift surveys from redshift space to real 
space, and in the process predict the linear cosmological 
velocity field and the Local Group dipole. In doing so it is 
useful to have some idea of what uncertainties arise in the 
process. In this Section we calculate some of the properties 
of the reconstructed density and velocity fields. Throughout 
we shall assume that the survey is spherically symmetric, 
which is nearly the case for the IRAS Point Source Redshift 
Survey (Saunders et al 1996). 

It is useful to disentangle the various uncertainties in 
analysing redshift surveys. In the following sections we treat 
independently the effects of: a bias ed and uncertain input 
distortion parameter, /3 (Section |4.l[ ); a finite survey volume 
on the reconstructed density field (Section 4.2); the proper- 



ties of the velocity field in the CMB frame for finite survey 
volu mes (Section 4.3); a finite survey on the dipole (Sec- 
tion 4.4); the properti es o f the velocity field in the Local 



4.1 



and 



4.2 



Group frame (Sec tion 4.5); shot-noise (Section it). Apart 
from Sections 



we shall ignore the effects of the 
redshift distortion. While this over-simplifies the case, it is 
useful to understand each of the effects in isolation. 



3.6 Filtering redshift surveys 

A number of methods for reconstructing the density and ve- 
locity fields do not use a smoothed galaxy distribution, but 
rather treat the galaxies in the catalog as test particles, in- 
verse weighted by the selection function. The question then 
arises, should one smooth or not? Clearly information is lost 
whenever one smoothes. However, there are a number of rea- 
sons why one can smooth in reconstruction. Firstly we are 
applying linear theory to the density field, so the scales we 
are considering must be large. Secondly, as well as nonlinear 
features in real space, there are nonlinear caustics and "fin- 
gers of god" in redshift space that we must remove before 
applying the present formalism. An alternative to smooth- 
ing "fingers of god" is to collapse them back onto the cluster 
they originate from. The small distances between galaxies in 



4.1 The wrong (3. 

Using the incorrect distortion parameter, (3, will lead to a 
systematic offset in the recovered fields. If the true distortion 
is f3, and we use a different value f3' in the reconstruction, a 
residual term is left in the reconstructed density field. This 
can be characterised by finding the product S'~ 1 S, where 
S'^ 1 is the inverse operator with the incorrect distortion 
parameter. This residual is 



S" 



■0 



j(d 2 r + 2r- 1 d r )V- 



(74) 



l+(3' 

for a constant offset distortion. Another possibility is that 
there is a large uncertainty on the distortion parameter (as 
is the case for the present generation of redshift surveys) 
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and a stochastic term is added to the reconstructed fields. 
If we now take j3' to be an unbiased but uncertain estimate 
of the true distortion parameter, the induced uncertainty on 
the density field is 



a si = 



on 



(1 + /3) 2 



(dr +2r- 1 d r )X7- 2 5. 



(75) 



Hence the uncertainty due to errors in the distortion param- 
eter lead to a density-dependent uncertainty in the recov- 
ered fields. 



4.2 Incompleteness in reconstructed density fields 

The reconstruction of the cosmological density field is always 
incomplete, since the distortion is nonlocal and we are trying 
to reconstruct with only a finite volume. Hence uncertainties 
arise in our estimates of the velocity and shear fields which 
cause the distortion. We can quantify this uncertainty in 
a model-dependent fashion by calculating the effects of a 
finite survey and finding the rms contribution from external 
structure not included in the survey volume. 

The effect of a finite survey volume is to alter the true 
inverse redshift operator, S , to an effective one, S c ff , using 
only the information within the surveyed region. Using the 
definition of the inverse redshift operator (equation EOJ) we 
see that this is equivalent to the operator 



sj = (v^v^v^ 2 

where V 



(76) 



fl is the inverse Laplacian defined over a finite vol- 
ume. It is useful to split the inverse Laplacian defined over 
all space, V~ 2 , into a term defined within the survey vol- 
ume, V~^, and a term defined outside the survey volume, 

V- 2 = V-^ + V;^. (77) 



By doing so we can find the difference between the recon- 
structed density field and the true real space density field 



terms of V 



equivalent to V 

vr ff 2 = v 



>« 2 

<Ri 



<R- 



Since our effective inverse Laplacian is 
we can write 

(78) 



Let us define 
A = [SjS]S-S 



(79) 



as the residual field after reconstruction. We find that this 
field can be expressed as 



-[s _1 v 2 v;^ 



s] s. 



(80) 



Expanding the density field in spherical harmonics and ap- 
plying these operators we find 



A(r) = A - 



P 



7T 1+/3 



dkk 2 5i m {k)- 



\{i+l)F t {k,R) 
21 + 1 {kr) 2 - e 




40 60 
r (h _1 Mpc) 
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Figure 3. Cosmic variance in reconstructed density field due to 
external structure for a volume— limited redshift survey with ra- 
dius R = 100/i -1 Mpc. The upper line is the total uncertainty 
from all modes. The line decreasing with radius is the 1 = 1 
dipole contribution. This diverges at small radius. The flat light 
grey line is the 1 = 2 quadrupole contribution. The increasing line 
is the total contribution from all modes i > 3. 



and diverges at the origin, suggesting that this should be 
removed in reconstructions by working in the Local Group 
frame (we shall return to this point when reconstruction 
velocity fields), while the 1 = 2 quadrupole residual is inde- 
pendent of radius. 

The variance of these terms can be found by squaring 
and ensemble averaging using the result that the ensemble 
average of the 8i m (k)'s is (Heavens & Taylor 1997) 



{Sem(k)S* e , m ,{k')) = P(k)k~ 2 S D (k - k')df e ,d. 



K 

mm' 



(83) 



where 8d{x) is the Dirac delta function and Sfj 
necker delta. The cosmological variance introduced by in- 
completeness is given by 



<A 2 (r)> = 



P 



1 + P 



D k 2 dk sr^ne+ir 



F e (k,R) 
(kr) 2 ~ e 



Y tm {r)(81) 



where 

F e (k,R)= / ^-(k 

J R 



r>?\i+ P - P ^mm.^) 



k 2 r'' 



(kr'Y 



For (3 = we have the trivial solution that the residual field 
is zero, since no reconstruction is necessary. We also note 
that the monopole mode is zero from Newtons' First Theo- 
rem. The £ = 1 dipole term increases with decreasing radius 



Figure g shows the radial dependence of the rms of 
this residual field for (3 = 1 and a maximum radius of 
100/i _1 Mpc. We again use the Peacock-Dodds form for the 
power spectrum. The upper dark line is the total rms resid- 
ual field. The overall impression is that for reconstructing 
density fields the uncertainty is low, except at the center and 
edges of the survey. Linear theory then suggests that den- 
sity fields can be recovered to high precision, ignoring the 
effects of a selection function and shot noise. To understand 
where this uncertainty comes from we have also plotted the 
rms residuals for the dipole and quadrupole. As expected 
the dipole increases at the origin, while the quadrupole is 
uniform across the reconstructed survey. The higher modes, 
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l > 3 contribute to a term increasing with radius, but which 
does not diverge at the survey boundary. Interestingly the 
dipole, quadrupole and the contribution from the rest of 
the multipoles are all equal at about half the survey ra- 
dius, independent of the survey radius. Thus a nice way to 
characterise the uncertainty is in terms of the quadrupole 
contribution, which is independent of radii, and a third of 
minimum uncertainty. 



4.3 Incompleteness in the velocity field 

The incompleteness that arises when reconstructing the den- 
sity field is not the only worry when calculating cosmic fields 
from a finite region. In the rest of this Section we consider 
the effect of a finite volume on estimates of the velocity 
field in the CMB and Local Group frames. To simplify the 
analysis we shall assume no distortion, and only consider the 
uncertainty due to a finite survey volume. We shall also con- 
sider the effects of shot noise on the reconstructed velocity 
field. We begin by expanding the velocity field in spherical 
harmonics. 

The Newtonian Green's function in equation (Q) can be 
expanded in spherical harmonics 



(85) 



and is related to the velocity field by V(l/r) = —r/r 2 . The 
radial function Ki(r,r') is given by equation (^) 

The gradient operator can be conveniently decomposed 
into radial and transverse terms; 



V[KiYi m (r)] = [d r Kt]Yi m {r)-i^l + l))^- 1 !^*- )(86) 
where we have defined the vector spherical harmonics 



rYtr 



1 tm 



:(r x L)Y tm . (87) 



and L = — ir x V is the classical angular momentum oper- 
ator. With these definitions we find that the radial velocity 
component is given by 



v r (r) = 

where 
Ui m (k, r 



21+1 



dkk Stm(k)Uem(k,r) 



dr'r' 2 [d r K e (r,r')]j e (kr')Y tm (r). (89) 



The transverse velocity components are given by 
Vt ( r ) = J dkk 2 5 tm (k)W em (k,r) 



(90) 



where the kernel is a transverse vector field 
W lm (k,r) 



(2t + l)r 



dr'r' 2 K £ (r,r')jt(kr')Y em (r). (91) 



Using equation (g3|) for the expectation value of two 
harmonic modes we find the variance of the velocity field in 
terms of the power spectrum is 



(v 2 r (r))= / ^P(k)\U em (k,r)\ 2 , 



(92) 



where the effective window function is 

i 2 



\U lm {k,r)Y 



and 



51 (21 + 1) 



dr'r' [d r Kt(r,r')]ji(kr') 



R 



(93) 



°°1.2 



k 2 dk 



2tt 2 

with the window function 



P(k)\W im (k,r)Y 



\Wim{k,r)\' = 



^ (21+1) 



dr'r' 2 Ki(r, r')ji(kr') 



(94) 



.(95) 



The range of these radial integrals in the window func- 
tions determines the source of the uncertainties. One of the 
main uncertainties in reconstruction is the effect of the den- 
sity field external to the survey volume. In the next Section 
we study the effects of structure beyond the survey. These 
results can also be used to calculate the expected variance 
in the velocity field from shells or bulk motions. 



4-3.1 Incompleteness due to external structures 

The range of the radial integral in these equations depends 
on the region of space under consideration. If we are deter- 
mining the sampling variance due to structures external to 
the redshift survey, then the range is r > R, where R is the 
radius of the survey (here we shall only consider sharp edges 
to the survey. A more complete analysis will include weight- 
ing and the effects of an angular mask). In this case r < r' 
and so k,i — r /r' l+1 . Hence the radial integral in equations 
© and @ is 



dr'r' 2 Ke(r,r')ji(kr') = k 2 (kr) e 



i ( j t -i(kR) 

(kRy- 1 



(96) 



The radial derivative with respect to r, in equation (|93|) is 
trivial. 

The total sampling variance due to external structure 
is then 



£(r/R) 2 



-1) .2 

'31- 



(kR) 



(97) 



The summation is only non-zero from 1 = 1, since external 
structure cannot affect the monopole mode. 

Figure ^ shows the radial, transverse and total rms ve- 



locities, « rms = yj (v 2 ), for a survey with R = lO O/i 1 Mpc in 



the CMB and Local Group frames (see Section 4.5) frame. 
The upper set of three lines correspond to the radial, to- 
tal and transverse rms uncertainties in the velocity field 
and its components in the CMB rest frame. Interestingly, 
for a spherical survey the uncertainties obey the relation 
(v 2 (r)) < (v 2 v (r)} < (v 2 (r)). At the origin the uncertainty 
is about lOOkms -1 , which is the dipole uncertainty for a sur- 
vey of radius R — 100/t -1 Mpc (see Section 1.4). The main 
effect is a rise in uncertainty from the center of the survey to 
the edges. However, even on the survey boundary the error 
is finite. 

The radial and transverse velocity cross-correlation 
functions are 

(vr(r)vr(r')) = 
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Figure 4. Cosmic variance in reconstructed velocity field due 
to structure external to a redshift survey with radius R = 
100/i — iJVIpc. The upper set of three lines are the radial, total 
and transverse rms uncertainty in the velocity field and its com- 
ponents in the CMB rest frame. The lower set of lines are for the 
same survey, but calculated in the Local Group rest frame. 



2tt 
and 

(vt(r)vt(r')) 



(2£+l) \R 2 



1-1 

r -^) j 2 e _i(kR)Vi(n)(98) 



<lk 



^ £(i + l) ( rr' 



t-~> (21 + i) V R 2 



— jU{kR)V t i 



where Vtip) is the Legendre function and fi = r.r'. No cross 
terms arise in the case of an all-sky survey, as we have chosen 
orthogonal projections. 



4-3.2 Sampling variance from interior structure 

Over the range r < R the radial integrals in equations ( |93| 
and (B5|) can be again calculated; 



dr'r' 2 Ki(r,r')ji(kr') = 
] ji(kr) 



( 2t+l \ 
\ k 2 J' 



k 2 j jy ' k 2 (kRy- 1 

Differentiating with respect to r yields 



(100) 



(101) 



In the limit that R — > oo we recover the expression for the 
true peculiar velocity field (Regos & Szalay 1989) 



TV ' 



'r(r) = \ — ^ dkk 2 8i> m (k)k 1 j't{kr)Y lm (r) 



(102) 



Figure 5. Uncertainty in the reconstructed dipole due to cosmic 
variance from external structure (dark line) and shot noise due 
to the discr ete sampling of galaxies in the survey (light line; see 
Section (l.tj ). The upper curve is the total rms uncertainty due to 
cosmic variance and shot noise added in quadrature. The power 
spectrum is described in the text and the redshift survey selection 
function is modelled on the PSCz (Saunders et al 1996). 



where we a dash on the spherical Bessel function denotes 
differentiation with respect to the argument. The autocor- 
relation function of the radial velocities is 



f°° rik 

Mr)tv(r')> = J ^P(k) 



t 



(2i+l)j' e (kr)j' e (kr')Te(ri(W3) 



and their variance is 

dk 
2ff 



(vr(r))=j ^P(k)J2(2e+l)j' e 2 (kr). 



(104) 



4.4 Incompleteness in the cosmological dipole 

As r — > only the I — 1 (dipole) term survives and we 
recover the cosmic variance on the dipole 



= J ^P(k)jl{kR). 



(105) 



In Figure |B| we show how the dipole uncertainty changes 
as a function of survey radius, R, smoothed on a scale of 
5ft _1 Mpc. We have again assumed the linear power spec- 
trum suggested by Peacock & Dodds (1994). At zero radius 
the uncertainty tends towards the 1-d rms velocity. This 
fairly well matches the observed local group velocity if one 
assumes that our motion is a fair estimate of the 3-d rms 
velocity. In that case K rms .id = i> rm s,3d/\/3 w 350kms _1 . 
As the survey radius increases the dipole uncertainty de- 
creases. However, in a flux limited redshift survey the uncer- 
tainty from shot noise increases with radii (see Section 4.6), 



and eventually dominates over the sampling variance. For 
the PSCz we find that the uncertainties are equal at about 
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R = 150fe~ 1 Mpc , thereafter becoming shot-noise dominated 
(see Section 4.6). 

For a power-law spectrum, P(k) = Ak", the integral 
can be evaluated and we find 



v£(0)> = 



-R- 



n Y(n — 1) cos nn/2, 



(106) 



2" + 1 TT 

where the spectral slope is in the range — 1 < n < 1. 

The equations for the variance and correlations can be 
also be simplified for power-law spectra. In this case the 
Bessel integral obeys the scaling relations 



dkk"j(_ 1 (kR) = 



' dkk n i 2 (kR) r{e+{n - 3)/2)T({5 - n)/2) (107) 
dkk Jo (kR) r{e + (3 _ n)/2)v({n _ 1)/2) • (107) 

Hence the cosmic variance can be evaluated in terms of the 
dipole uncertainty; 



(Vr(r)) 



3<«, 2 (0)> 



\ " 



e 



^21 + 1 \R 



/ r \ 2(e-i) 



r(l+(n-3)/2)r((5~n)/2) 
r(^+(3-n)/2)r((n-l)/2)' 



(108) 



The transverse and total variances can be expressed simi- 
larly. Hence we find that the uncertainty on the velocities in 
the survey only scale as the dipole uncertainty for a power- 
law spectrum of fluctuations. In addition the uncertainty 
from the error in the dipole can simply be subtracted by 
removing the I = 1 mode for an all-sky redshift survey, cor- 
responding to moving to the Local Group rest frame. In the 
next section we find that this is generally true, for arbitrary 
power spectra. 



4.5 Transforming to the Local Group rest frame 

As we shall see it is more accurate to calculate velocities 
in the observer's Local Group rest frame, rather than the 
CMB frame used so far. Projecting the local dipole along 
and arbitrary radial direction we find that 



IV (0) = v(0).r = J— dkk 2 5 10 (k)k- 1 jo(kR) (109) 




Similarly the transverse dipole projected perpendicular to 
an arbitrary radial vector is 

/ i r°° 

#) = 1/r» / dkk 2 {Su{k)-5 10 (k)r i )k- 1 j (kR) (110) 
V 6tt 2 J q 

Similar results for the dipole have been found by Webster et 
al 1997. 

With these equations it is straightforward to show that 
subtracting the dipole term is equivalent to ignoring the 
monopole and dipole terms in the harmonic summation over 
£. This is because the local group velocity and other veloc- 
ities are only connected by a dipole term. In Figure ^ we 
plot the uncertainty in the velocity field in the Local Group 
rest frame. With the major source of uncertainty, the dipole 
mode, removed we find that the relative uncertainty in the 
velocity field is small compared with the magnitude of the 
flow fields we expect. 



4.6 Shot— noise contribution to the reconstructed 
redshift survey 

The shot-noise contribution to the density field can also be 
calculated from equation (^), taking into account that the 
autocorrelation function for Poisson noise is 



(S(r)6(r')) = ^-8 D (r-r'), 



(111) 



where <j>(r) is the survey selection function. We find that the 
shot-noise variance is (Taylor & Rowan-Robinson 1993) 



R '2 i / 

r dr 
(h(r')\r 2 — r' 2 \ 2 



(112) 



However, this diverges as r approaches r' , as infinite variance 
is produced by infinity close discrete particles. To avoid this 
problem in the reconstruction one smooths the density field. 
Applying a Gaussian filter here we find that the smoothed 
shot noise variance is 



where 



R '2 . / r 1 
r dr 



dn\G(r - r')\ 2 



+i 



G(r) = rC- 



rt{r/V2R s ) 



2 e" 



7T rR s 



(113) 



(114) 



is the Gaussian smoothed Greens function. 

The shot noise uncertainty on the velocity can also be 
decomposed into radial and tangential components along the 
line of sight. The radial term is 



sn (r)> 



R J2 



dr' 



dfi\r.G(r — r 



(115) 



and the tangential component can be found from the relation 



Figure q shows the shot noise contribution to the dipole 
for the example of the PSCz. We see that the shot-noise 
contribution rises just as the sampling variance contribution 
is falling (Figure |l|) . 

Figure 3 shows the radial, transverse and total rms ve- 
locity dispersion due to shot noise in the CMB and Local 
Group rest frames. The overall effect of moving to the Lo- 
cal Group rest frame is marginal, and only significant for 
velocities at small radii. Hence moving to the Local Group 
rest frame does not affect the shot-noise contribution to the 
uncertainties. 



5 CONCLUSIONS 

In this paper we have derived the inverse redshift space op- 
erator. By applying this operator to a redshift survey with 
the correct distortion parameter, found either from some 
other method, eg a study of the anisotropy of clustering in 
redshift space, or by comparing the reconstructed velocity 
field with the true one, a real space map of the density field 
can be recovered. From this one can reconstruct to linear 
order the real space potential and velocity fields. We have 
derived inverse operators for redshift surveys for observers in 
the CMB or Local Group rest frames, and we have derived 
operators which include the effects of a selection function. 
As a corollary to the calculation of an inverse operator 
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Figure 6. Shot noise in reconstructed density field due to in- 
ternal galaxy distribution in a redshift survey with radius R = 
100/i _1 Mpc. 



in the Local Group rest frame we have also shown how to 
calculate the observer's dipole directly from a redshift sur- 
vey, without having to fully reconstruct the density field, 
again in either the CMB or Local group rest frames. Simple 
tests on an ensemble of Gaussian random density fields have 
shown that both the reconstruction of the density fields, and 
the observers dipole are accurate. 

To simplify calculations in redshift space we have de- 
veloped the formalism using a spherical harmonic repre- 



sentation tor the fields and the distortion operators, this 



approach allows one lo easily esLimale the reconslmcLed 

field from a spherical harmonic decomposition of the redshift 
space density field. We have also found expressions that re- 
late the real space dipole to the redshift space density dipole. 

The spherical harmonic representation is also used to 
estimate the effects of a finite survey volume on the recon- 
struction. We find that the uncertainty on the dipole mode 
diverges at the origin, and suggest that this can be removed 
by working in the Local Group rest frame. The quadrupole 
uncertainty is a constant across the survey volume, while 
the sum of the remaining multipoles increases with radius. 

The velocity field can also be inferred from a recon- 
structed redshift survey and we have calculated the effects 
of a finite survey volume and shot noise on the reconstruc- 
tion of the velocity field, for the special case of no distortion. 
We have calculated the sampling variance expected for the 
velocity field reconstructed from a redshift survey and find 
that the velocity field relative to the observers motion can 
be calculated far more accurately than absolute velocities. 
Hence one should work in the Local Group rest frame when 
comparing velocities. Calculation of the shot-noise contri- 
bution shows that the shot-noise is fairly insensitive to the 
rest frame. 
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